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Abstract. We present the results of far infrared (FIR) spec- 
troscopy of the luminous blue compact galaxy (BCG) Haro 
1 1 (ESO 350-IG38) obtained with the Long Wavelength Spec- 
trometer (LWS) on the Infrared Space Observatory (ISO) in 
low resolution mode. This metal poor dwarf merger is an ex- 
tremely hot IRAS source. We discuss the balance between dust 
and line cooling in the photodissociated regions (PDR), in par- 
ticular the role of the [CII]A158/i line, and derive the basic 
properties of the PDR gas and estimates of the gas and dust 
masses. The mass ofthePDRs, 2t 2 10 8 M o, is comparable to 
that of the ionized gas and exceeds the observed upper limit of 
the HI mass. The gas/dust mass ratio is low, indicating that the 
galaxy contains little cold dust. The low metallicity, the intense 
radiation field and the low column density of Haro 1 1 results in 
an extremely high [CII]/CO flux ratio and probably also a very 
high A^(H2)/Lco conversion factor. Therefore CO is a poor 
indicator of the H2 mass in starburst dwarf galaxies. 

After a reanalysis we confirm the claimed correlation (Mal- 
hotra et al. |1997| ) between the [CII]158^/FIR flux ratio and 
the IRAS f6o/fioo dust temperature and reduce the scatter. We 
find that Haro 1 1 deviates from the relationship being brighter 
in [CII] than what would be expected, if the mechanism pro- 
posed by Malhotra et al. is dominant. As an alternative (or com- 
plementary) explanation we propose that the [CII]158/j/FIR 
versus f6o/fioo relationship is caused by an increasing optical 
depth with increasing IRAS temperature. The low metallicity 
of Haro 1 1 and its extreme starburst properties probably allows 
the medium to be thin despite its high feo/fioo ratio. This leaves 
room for a more optimistic view on the possibilities to detect 
massive starforming mergers at high redshifts, using the [CII] 
line. 
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1. Introduction 

Blue compact galaxies (BCGs) are characterized by their high 
star formation rate and low chemical abundances. The high ra- 
tio between blue luminosity and HI mass, Lb/Mhi, combined 
with a high relative HI mass fraction, indicates a high star for- 
mation efficiency and a high gas consumption rate. The appar- 
ent conflict between intense star formation and low chemical 
abundances can only be explained if the gas containing the re- 
cently produced elements is expelled from the galaxy, if fresh 
gas is accreted onto the galaxy or if the galaxy is young (e.g. 



Searle & Sarg ent |1 972^ Makarova & Karachentsev |1998| ; Pa- 
paderos et al. 1998). We still have rather vague ideas about 



what is triggering the intense global burst. Our previous obser- 



vations (Ostlin et al. 1999a 1999b) indicate that gas dynamical 



instabilities in connection with mergers could be the key mech- 
anism. Under such circumstances the heating processes of the 
ISM could become a complicated mixture of shocks caused by 
infalling clouds and winds from massive young stars and on the 
other hand a strong radiation field from the starburst. 

While optical and radio observations give us important in- 
formation about the conditions in the hot ionized and the cold 
neutral gas in star forming galaxies we often lack information 
about the the ISM in the transition zone, the photodissociation 
regions (PDRs). The gas in PDRs contains neutral atoms and 
ions of low ionization potential {\ < 13.6 eV), e.g. C + . With 
the different possible heating mechanisms involved one would 
expect that the geometry and spatial extension of the ISM and 
the PDRs would be quite different in galaxies like Haro 1 1 and 
those with quiet star formation activity or with active nuclei, 
resulting in radically different spectral signatures. 

Spectral diagnostics are valuable tools in dissecting the 
balance between different heating and cooling processes and 
the physical conditions in the gas. The first observations of 



[OI]63^in M82 was presented by Watson et al. (1972b). Craw 



ford et al. (1985) for the first time showed that important in- 



formation about the PDR physics in external galaxies can be 
obtained from the collisionally excited [CII]A158/i. Later, sev- 
eral studies of star forming regions in the Milky Way as well 
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as other starburst galaxies and ultraluminous infrared galax- 
ies ('ULIGs', possibly hiding active nuclei) have demonstrated 
the power of more detailed far-IR spectroscopy (e.g. Colbert 



Optical spectroscopy was obtained with ESOs 1.5-m tele- 
scope in 1983, using the Image Dissector Scanner and an aper- 
ture of 4"x4" arcseconds 2 (Bergvall & Olofsson 1986). Addi- 



et al. 1999). An advantage with observations in the IR region 



is that we diminish the influence of dust extinction, allowing a 
more direct comparison between 'naked' starbursts in low mass 
galaxies and hidden starbursts in ULIGs. 

The fact that [CII]A158/i emission constitutes the peak of 
the flux distribution f v in starforming galaxies has led to a dis- 
cussion about the possibility of using [CII] to detect young 
galaxies at high redshifts (e.g. Petrosian et al. 1969 and the de- 
tailed discussion by Stark 1997). Despite some emerging pes- 
simism about this possibility (Gerin & Phillips 1998) we will 
argue that this still is a viable idea. 

As part of our multifrequency study of three of the most 
luminous BCGs in the southern hemisphere this report will 
focus on observations u sing t he Long Wavelength Spectrom- 
eter (LWS) (Clegg et al. 1996 ) on the Infrared Space Observa- 
tory (ISO) (Kessler et al. |l996| ) of the luminous BCG Haro 1 1 . 
Only a few emission lines in the LWS window are detectable in 
most nearby BCGs in a reasonably long integration time. In our 
case we chose to obtain clear detection of a few of these lines, 
with a sufficiently high S/N to allow comparisons with models 
and other galaxies. These lines are the [OI]63/i, [OIII]88/i and 
[CII] 158/1 lines, which constitute important coolants in the far- 
IR regions. Among these lines, [CII] is of particular importance 
for reasons just mentioned and discussed below. In addition to 
these three lines, we obtained an upper limit on the flux in the 
[OI] 145 fi line in one of the adjoining detectors. 



2. Observations and reductions 

The observations were carried out during the routine observa- 
tion phase of ISO in the low-resolution LWS02 observation 
mode with a resolution of 0.6 p. We used a scan width of 7 
resolution elements, 8 samples per element and fast scanning. 
Each spectrum was scanned 6-111 times with a nominal 0.50 s 
integration time giving in total an effective integration time of 
2820 s. 

The data presented in this paper have been taken from 
pipeline 6.1. All the glitches in each individual scan have been 
removed by using the Interactive Analysis Software Package. 
About 15% of the total scans were rejected due to heavy cos- 
mic rays contamination. In addition to the requested data the 
LWS spectrometer provided data on all the other detectors. 
These data have been also analyzed and an upper limit on the 
fine structure line of [OI]A145/i has been derived. The result- 
ing data were resampled and coadded. The data are shown in 
figure 1. 

Exept for the LWS data we also present some additional 
previously unpublished data in Table 1 . The CO observations 
were carried out during June 9-12 1988 with the SEST tele- 
scope equipped with coole d Schottky receivers at 115 GHz 
(Jorsater & Bergvall 1999 ). Six blue compact galaxies were 
observed and none was detected. 



tional data were obtained with the IDS at the ESO 3.6-m tele- 
scope in 1986. Standard reduction techniques were used. More 



information is found in Bergvall & Osttlin (1999a) 



3. Results 

3.1. General properties of the target galaxy 

Haro 11 (ESO 338-IG04) is a galaxy with extreme starburst 
properties, summarized in Table 1 . In the following we will as- 
sume a Hubble constant of Ho = 65 km s _1 Mpc _1 . Most of the 
data are from our pre vious observation s in the UV-rad io region 
(Bergvall & Olo fsson |1986| ; Ostlin et al |l999a| , |l999bfc Bergvall 
& Ostlin 1999a and other unpublished data). Our observations 
show three conspicious condensations in the ce ntre e mbedded 
in a regul ar Ho halo. HST images (Malkan et al. 1998 ; Ostlin et 
al 1999b ) reveal a merger morphology. Faint whisps and indi- 
cations of shell structures characteristic of galactic disk progen- 
itors are seen in the outer regions. Our studies of the kinematics 
of the ionized gas indicate that the starburst is triggered by coa- 
lescing clouds involved in the merger process. The Ha velocity 
field is requires a mass of the main body of >10 9 Jv[ Fj (Ostlin 
etal. 



1999a; 1999b) 



The oxygen abundance in the central region was calculated 
from the optical spectra, using the standard method. We as- 
sumed a homogeneous single metallicity model with two tem- 
perature zones. The electron temperature of the ++ zone was 
derived from the [OIII] A4363/([OIII] AA(4959+5007) ratio and 
that of the + zone from the semi-empirical relationship de- 



rived by Vila-Costas and Edmunds (1993 ). The density was de 



rived from the [SII]A6717/6731 line ratio and was found to be 
low. With this approach we obtain a metallicity close to 10% 
solar. Such a value is unusually low for a ga laxy o f this lu- 
minosity and deviates strongly (Bergvall et al. 1998 ) from the 
metallicity-lumi nosity relationship of normal dwarf galaxies 
(Skillman et al. 1989 ). Part of the optical emission may have 
a shock origin and we made an estimate how large the con- 
tribution from shocks may be by mi xing p redicted line spec- 
tra from shocked regions (Raymond |1979| ) with those of HII 
regions from Stasihska (1982). The best fits generally gave a 
contribution from the shock component to H/3 of < 5%. This 
is in most cases not enough to change the result of the calcula- 
tions of the metallicities drastically. 

A surprising fact is that no HI has been detected in 
the galaxy (Bergvall et al. 1999b). With an upper limit of 
Mhi =10 8 Mq and consequently MHi^B,totai <0.01, it 
seems to be remarkably devoid of neutral hydrogen for its star 
formation rate. One possible explanation is that the neutral gas 
is ionized by the central starburst. Another possibility is that 
much of the gas is in molecular form. The ISO observations 
discussed here, which relate to the transition zone between the 
ionized and the cold molecular medium, will help to derive a 
coherent picture of the different phases of the gas. 
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Fig. 1. The ISO LWS emission lines of Haro 11.3 detections are presented. An upper limit of the flux of the [OI]A145/i line was 
derived from the region at the expected position of the line, marked with an arrow. The bars indicate mean errors per pixel. 



3.2. The dust properties 

From the 60 and 100 p IRAS fluxes and assuming a FIR grain 
opacity of k = 2.5 10 3 A -1 cm 2 g~ 1 (Hildebrand |l983| ), where A 
is the wavelength in microns, we calculate a dust temperature 
of Tdust = 49K. We can then obtain a rough estimate of the 



mass of the warm dust from (Devereux & Young 1990) 



M 



dust 



4.58 / 100 r 2 (e 144 / T - 



DM, 



fioo is in Jy and r is the distance in Mpc. We obtain 
M du st ~ 4 10 6 M Q . 

It is interesting to compare gas-to-dust mass ratio with that 
of other gas rich galaxies. In Sect. 3.4.3 we find that the total 
hydrogen gas mass, i.e. HI, HII and H2, is approximately 10° 



M P) . We thus obtain M gas IM 



l dust 



< 200. In a study of normal 



spiral galaxies, Devereux & Young ( 1990| ), taking into account 
only the inner part of the HI disk containing the warm dust, 
found that M ga slMdust ~ 1100. There is thus a difference a 
factor of 5-6 difference between Haro 1 1 and normal spirals. If 
the dust content correlates with metallicity, one would expect 
that M-gaslMdust would increase with decreasing metallicity 
and indeed this i s the case for dlrrs as was shown by Lisenfeld 
& Ferrara (1998). If we apply this metallicity correction in the 



comparison, it would lead to an even larger difference between 
Haro 11 and the spirals in the Mg as IMdust ratio. The high 
value for spiral galaxies has been shown (e.g. Chini et al. 1995) 
to be due to the fact that a large amount of the dust is of the 
cold, 'cirrus' type, radiating at longer wavelengths and thus not 
contributing to the mass estimates based on the IRAS data. The 
fact that Mgasl Mdust of Haro 1 1 is so low (in fact close to the 
value found for our own Galaxy, including the total dust mass, 
i.e. MgasIM dust, total ~ 100-150) indicates that the galaxy 
contains very little dust at low temperature. 



3.3. The LWS spectra 

Figure 1 shows the spectral profiles of the three detected lines 
and the region around [OI] 145/i, for which we derived an upper 
limit. Table 2 summarizes the spectral data. From our optical 
spectra we have derived an extinction coefficient in the bright- 
est region of the galaxy of Ay ~ 0.8 mag. Thus the corrections 
for galactic extinction in the infrared are considerably less than 
the observational errors in the observed emission lines (Adams 



et al. 1988) 
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Table 1. General properties of Haro 1 1 



Ref./Note 



01950 


00 h 34 m 25 s 7 




<5l950 


-33°49'49" 




m fl 


14.3 


1/a 


vo 


6175±10kms _1 


1/b 


r 


92Mpc 


3 


M s 


-20.5 


a 


a MB =26.5 


9kpc 


2/c 


B-V 


0.5 


1/a 


12+[0/H] 


7.9 


2 


^Lya. 


9 10 35 W 


3 




3.2 1 35 W 


2 


Lfir 


7.4 10 10 L 




feo/fioo 


1.48 




Lfir/Lb 


5.7 




Loo(i-o) 


< 1.0 1 29 w 


4 


logCLco/LF/ij) 


<-8.5 






< 10 s M Q 


d 




< 10 s M Q 


4 




< 0.01 





a) Corrected for galactic extinction (Burstein & Heiles 1982) 

b) Radial velocity, corrected to the local group according to 
vo=v+300 sinZ cosfe, where v is the heliocentric velocity and / and b 
are the galactic coordinates 

c) r: distance; a: major axis diameter 

d) Us ing a conversion factor of 3 10 cm" ~ 2 (Kkms -1 ) _1 (Sanders 
et al. 



1991) 



Refere nces: 1) Bergvall & Olof sson |1986| , 2) Bergvall & Ostlin 
19994 3) Bergvall & Olofsson, |l999c[ , 4) Jorsater & Bergvall, [l999| 



Table 2. Observed LWS data. The table shows the measured 
line flux, its estimated mean error and equivalent width (W) 
and the adjacent continuum level flux density. 



Line id. 


Flux 


W 


Continuum flux density 




(10 _16 Wm- 2 ) 


to 


(10- 16 Wm- 2 /i" 1 ) 


[OI]63aj 


9.2 ± 1.8 


.22 


42 


[OIII]88/i 


25 ± 0.8 


.84 


33 


[OI] 145 ^ 


< 0.3 


< 0.2 


1.7 


[CII]158m 


3.9 ±0.1 


2.0 


2.1 



3.4. Properties of the PDR gas 

One of the purposes of this study is to investigate how much 
the PDR gas contributes to the total mass in this type of ob- 
jects, and in particular Haro 1 1 due to its remarkably low HI 
content. It is known that the relative volume of the PDR that 
occupies the weakly ionized (only low ionization stages) gas 
may be quite substantial. Regions having a total extinction of 
Ay < 10 are partly photodisso ciated (see e.g. Crawford et al. 



about the PDR physics). This is the situation for most of the 
neutral gas in the Galaxy. Estimates also show that < 40% of 
the total gas mass in starburst nuclei are PDRs, i.e. the pho- 
todissociated gas is the dominant phase in these extreme en- 



vironments (Stacey et al. 1991 ). The most widely used model 



for calculating the energy balance and obtain predictions about 
emissivities of the gas and dust was presented by Tielens & 



Hollenbach (1985). Later upgrades and applications have been 



discussed by Ho llenba ch et al. (1991), Womreet al. (1990) and 
Kaufman et al. ( |l999| ). 

The Kaufman et al. standard model is based on an oxy- 
gen abundance of 40% solar, close to the local Galactic value 
and a factor of 4-5 higher than in Haro 1 1 . From our derived 
carbon abundance in ESO 338-IG04 (Tol 1924-416), a galaxy 
quite similar to Haro 1 1, we may assume that N(C)/N(0) « 0.2 
(Bergvall 1985). Similar values have been obtained for other 



metal poor galaxies (Garnett et al. 1995). Thus, oxygen and 
carbon are underabundant with factors 4-8 with respect to the 
standard setup in Kaufmans et al. models (their models have 
N(C)/N(0) sa 0.5). The predicted temperature of the surface 
gas layers obtained from their basic diagnostic diagrams is 
therefore probably somewhat overestimated while still the line 



intensities are reasonably correct (Wolfire et al. 1990). How- 
ever, Kaufman et al. also discuss a case of lower metallicity, 
about 4% solar and present diagrams showing how [CII], FIR 
and CO luminosities relate to density, n, incident far UV ra- 
diation intensity, Go, and the two metallicities w 4 and 40% 
solar. In the following we use these models to estimate some 
fundamental parameters of the PDR regions. 

3.4.1. Optical thickness 

To enable a comparison of the line intensities with the mod- 
els of Kaufman et al., we first estimate the contribution from 
HII regions to the line intensities and also the optical depth 
of the lines. The optical emission lines (Bergvall & Olofsson 



1986 and newer data) have been used together with recent pho- 



toionization models (Stasinska 1990, and updates) to estimate 



1985| or Tielens & Hollenbach |1985| for a general discussion 



the contribution to the far-IR lines from the HII regions. We 
used our observed line intensities in [OI], [Oil], [OIII], [Nil], 
[SII], Hel and Hell, corrected for extinction, underlying ab- 
sorption and shock contribution to select the best fitting HII 
region models from Stasinskas library. As a result we find that 
the measured intensity of the [OHI]A88/i line is fully consis- 
tent with photoionization by a cluster of stars with T e f f of the 
order of 35000-40000 K. For the other 3 lines observed with 
the LWS, the contribution from HII regions to the measured 
intensity is strongest in [CII] but is still < 20%. Thus, when 
discussing the [CII] line flux and flux ratios, we can directly 
compare to the predictions from Kaufman et al.. 

In a study of two star forming complexes in the Large 
Magellanic Cloud, Israel et al. ( 1996| ) concluded that these 
regions were optically thin at a minimum column density of 
3 10 21 cm~ 2 . Since these cloud aggregates probably have sim- 
ilar properties as Haro 1 1 we can conclude that the star forming 
regions in Haro 1 1 on a local scale probably are optically thin 
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in 158/i. The optical extinction derived from the Ha/H/3-ratio 
in Haro 1 1 is only Ay =0.8 so we don't expect that the [CII] 
line is optically thick even on a global scale. But the derived 
extinction depends on the spatial distribution of the HII regions 
so to assure that we don't have conditions that allow hidden 
sources in the central area we will estimate the optical depth 
in [CII], Ten- This calculation will also be used to make es- 
timates of optical depths of the other galaxies involved in the 
discussion. 

In the high-temperature (T^92K, the excitation potential 
for [CII]), high-density limit (n^>4 10 3 cm~ 3 , the limit for col- 
lisional deexcita tion b y H and H2), t\cii] can be expressed as 



(Crawford et al. 1985) 



T[CU] 



2.3 1(T 16 



N c+ 
Tc+cr v 



where Nc+ is the column density of C + in cm~ 2 , Tc+ is 
the temperature of the gas in K and <j v is the velocity dispersion 
in km s _1 . If we as a first rough approximation assume that 
the global distribution of the gas is close to a homogeneous 
exponential disk we can derive the central column density in 



atoms cm 



N„ 



1.3 10 



14 



Mggs 



where A4 gas is the mass of the gas in solar masses and 
h gas is the scale length in kpc. As a first approximation we will 
assume h c .+=h# a . From our Ha images we derive h^a ~ 2" 
= 1 kpc. As an estimate of the gas mass we will use our upper 
observational limit of Mhi, i-e. around 10 8 Mq- With this we 
obtain a central column density in HI of Nhi ~ 2 10 21 atoms 
cm -2 . Later we will find it consistent to assume that most of the 
HI actually is associated with the PDR regions. Part of the C + 
region in the PDR also contains molecular gas. In sect. 3.4.2 we 
will show that this may increase the estimated column density 
with a factor of 2. 

The C + column density can now be derived from 



N, 



c+ 



1.3 10 14 X(C" 1 



M 



HI 



2-Kh 2 Ha 



where X(C + ) is the relative abundance of carbon ions to 
hydrogen atoms. With these approximation we can write the 
optical depth as 



r [C ii] < omx(c+)- r 



M 



HI 



' T c+ a v 2-Kh Ha 

Since most of the carbon in the [CII] emitting re- 
gion is in ionized form we will assume X(C + )~X(Ctotai)~ 
j^lgn X(0)« 1.2 10~ 5 . In starburst galaxies Tc+ is typically a 
few hundred K. We will later iteratively derive Tc+ and obtain 
a value below 400K. From our Fabry-Perot observations in Ha 
of the central region we estimate a v ~ 60 km s _1 . We then ob- 
tain t\cii\ < 0.001 . In order to derive a more realistic estimate, 
we have to take into account that the gas has a lumpy, fractal 
distribution. But it is sufficient to note that, assuming a volume 



filling factor a of a few per mille (Crawford et al. 1985| ), the 
surface density would increase with ~ a" 2 / 3 , i.e. a few tens. 
Thus it still would be safe to assume that the medium is thin in 
[CII]. The optically thick limit at a chemical abundance similar 
to Haro 1 1 would be at a mean column density of Nhi ~ 10 23 
atoms cm~ 2 and with solar abundances Nhi ~ 10 22 atoms 
cm -2 , still with the assumption that the filling factor is a few 
per mille. 

3.4.2. The mass of the PDR gas 

From the diagnostic diagrams by Kaufman et al. we can ob- 
tain the density and intensity of the incident far-UV flux Go 
in the PDR region. We use two diagrams for this purpose, one 
utilizing the [OI]63/i/[CII] 158/i intensity ratio and one utiliz- 
ing the ([OI]63/i+[CII]158/i)/FIR ratio. As mentioned above, 
before we use the predictions we need to correct for the differ- 
ence in metallicity between Haro 1 1 and the models, a factor 
of 3-4. We can estimate what effect this may have by compar- 
ing the predicted fluxes of these lines in an HII region at the 
different metallicities. The column densities and thus the vol- 
ume emissivities of the low ionization lines are lower in HII 
regions than in a typical PDR region but the flux ratios should 
vary in approximately the same way. From Stasihskas models 
we find that the difference in flux ratios, assuming a constant 
N(C)/N(0), is about 2 %. N(C)/N(0) probably differ with a 
factor of about 2 however, and so we should make an appro- 
priate correction for this which means a factor of 1.5-2 in the 
[OI]63/i/[CII] 158/i ratio. Kaufman et al. show in later diagrams 
that the metallicity has a minor effect on the [CII] 158/i/FIR ra- 
tio in the region of the parameter space where we are residing. 

We now feel confident to compare the slightly corrected 
line intensities with the predictions from the PDR models. Us- 
ing the [OI]63/i, [CII] 158/i line fluxes and the FIR flux, Kauf- 
mans et al. diagnostic diagrams show that Go is 1-2 10 3 in units 
of the Galactic far-UV field (1.6 10~ 6 W m" 2 ). The obtained 
value for the density is n « 2 10 3 cm~ 3 and thus Go/n 1-2 
cm 3 . These are typical values for normal starforming regions 
and galaxies (Tielens & Hollenbach, |l985[ Wolfire et al. |1990; 



Carral et al. 1994; Fischer et al. 1996). The surface temperature 
is T ~ 400K. The data are summarized in Table 3. With these 
data at hand we can obtain the predicted [OI] 145/i/[OI]63/x ra- 
tio which is approximately 0.06. This is within the errors of 
the measured line ratio and thus supports the reliability of the 
model. 

From the measured [CII] flux we may now calculate the 
mass of the warm atomic gas according to (e.g. Wolfire et al. 



199C) 



Mpdr = 5.96 10 15 



r 2 f[cii] m H 



Mr 



K{CII) X{C+, 

where r is the distance to the galaxy in Mpc, £\cn] is the 
[CII] flux in W m -2 , m# is the mass of the hydrogen atom in 
kg, and A(CII) is the cooling rate per atom in the 158/i line in 
W. In the high-temperature and high-density limit A(C II) wil l 
be nearly independent of the temperature (Wolfire et al. 1990). 
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A(CII) can therefore be rather safely estimated and we obtain 
A(CII) ps 1.3 10~ 26 W atom" 1 . In the surface layers of the 
PDR carbon is mainly in the low ionization monoatomic state 
so X(C+)ps X(C). We then derive M PD r~ 2± 2 10 8 M Q . 



3.4.3. The mass fractions of cold, warm and ionized gas 

The value we have derived for the mass of the PDR:s is larger 
than or comparable to our upper limit of the HI mass, i.e. 
10 s A4q, indicating that most of the HI gas is located in the 
PDR regions if these are dominated by the atomic gas com- 
ponent. The total molecular gas mass in normal star forming 
regions is normally 5-10 times larger than in the PDR regions. 
Data from other starburst galaxies however indicate that con- 



Table 3. Global data, based on the LWS observations 



siderably lower values are more likely (e.g. Wolfire et al. 1990; 
Stacey et al. 1991). This is probably a metallicity effect. In 
galaxies with low metallicities the UV radiation can penetrate 
deeper into the molecular clouds, increasing the PDR zone rela- 
tive to the CO core. An upper limit is therefore probably around 
10 9 M.q in our case. This mass estimate may be compared 
to the upper limit of the H2 mass, derived from the CO (1- 
0) observations which is < 10 8 Mq. The calculation of this 
last value is however rather uncertain due to the uncertainty in 
the conv ersion factor betwe en the CO flux and Ai h 2 (Maloney 
& Black |198S| , Taylor et al |1998[ ). As just mentioned, the low 
metallicity also allows a larger proportion of H2 to be hidden in 
the CO dissociated zone. Therefore, in addition to the problem 
caused by the uncertainty in the conversion factor between CO 
and H2, the CO fluxes are poor indicators of the H2 masses in 
low metallicity galaxies. 

Using the available photometry of the halo and the burst, 
our spectral evolutionary models predict a total stellar mass of 
1.5 10 10 A4q , while the observed rotation results in a mass 
of ~ 2 10 9 M.(?) , assuming dynamical equilibrium (Ostlin et 



al. 1999a ; 1999t| ). The discrepancy between the photometric 
and dynamical mass estimates can be resolved if the galaxy is 
dominated by velocity dispersion, which is reasonable if the 
width of the Ha line reflects potential motions, or if the galaxy 
is not in dynamical equilibrium. 

From the global Ha luminosity and the mean density of the 
gas, derived from our spectra, the mass of the ionized gas is 
estimated to be - 10 8 - 10 9 M G (Ostlin et al. |1999a|jl999b[ ). It 
is a remarkable fact that the HI mass is so low and also lower 
than or comparable to the estimated mass of ionized or molecu- 
lar gas. A possible explanation may be that HI to a large extent 
has become molecular and/or ionized in the merging process. If 
this is correct, then the gas/total mass estimates may have been 
severely underestimated in many dwarf starburst cases. This 
will complicate the efforts to understand dwarf galaxy evolu- 
tion in general. A summary of the mass estimates is found in 
Table 4. 



3.5. Thermal balance and the [CIIJ/FIR ratio 

The important coolant [CII]A158/i originates from the surface 
layers of PDRs where the C + ions are excited by photoelec- 
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Fig. 2. The IRAS f2s/feo versus the f2s/fioo micron flux ratios 
for Haro 1 1 (star) and the comparison sample (filled dots) 



trons ejected from grains heated by the UV radiation (Watson 
1972a| ). PAHs are believed to be the dominating agent of this 



process. Stacey et al. (1991) studied a mixed sample of spiral 
galaxies, starburst galaxies, giant molecular clouds and galac- 
tic star forming regions and noticed that while there is a linear 
relation between the energy density of the UV field and the 
cooling rate of the dust, the efficiency with which the [CII] 
line is excited seems to decrease with the intensity of the UV 
field. One should have in mind however that this is not a homo- 



geneous sample. Malhotra et al. ( 1997 ) investigated the same 
relationship between the [CII] cooling and dust cooling for dif- 
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ferent types of galaxies, using the [CII]/FIR luminosity ratio. 
An anticorrelation was found between the [CII]/FIR ratio and 
dust temperature, as measured by the IRAS f6o/fioo ratio. Af- 
ter a discussion of several possible explanations, Malhotra et al. 
found it most probable that this is a consequence of a decreas- 
ing efficiency of ejection of photoelectrons from the grains, re- 
ducing the [CII] intensity. Theoretically it can be understood 
if the ratio of UV flux to gas density, and thus the dust tem- 
perature, becomes high because the dust grains would be more 
positively charged, forming a stronger Coloumb barrier against 
photoejection. 

One should be aware, both as regards the results from Mal- 
hotra et al. and in the following discussion, that the sample 
for which useful data are available is quite heterogeneous. The 
galaxy types involved range from normal galaxies to Seyferts 
for which the relative importance of different heating sources 
(stars and AGNs) and different chemical abundances are not 
satisfactorily controlled. Moreover, in the first approximation, 
two dust components, warm and cold, having different spatial 
distribution, can be separated. This increases the complexity in 
the interpretation of the global data and the 60/100 index as 
well as the Lfj^/Lb index cannot be regarded as clearcut trac- 
ers of the star formation activity. As a consequence we expect 
that the scatter in the discovered trends will remain quite large. 

We will now compare the data for Haro 1 1 with a data set 
similar to that used by Malhotra et al. but including more de- 
tected galaxies at relatively high f6o/fioo- These data were ob- 



et al. 1994; Taniguchi et al. 1990), destroying the smallest dust 
grains (e.g. Desert |1990| ). 



tained from Stacey et al. (1991 ), Luhman et al. (1998), Malho 



tra et al. (1997) and Lord at al. (1996) 



For the most nearby galaxies, having velocities < 500 
kms -1 , we have scanned the literature to find the most re- 
cent distance determinations derived from photometric dis- 
tance indicators. This concerns the galaxies M82, NGC 3109 
(Kennicutt et al. Il998|), NGC 4414 (Turner et al. Il998|), 



NG C 47 36 (Garman & Young |1986[), NGC 5128 (Tonry et 
al. |1995| ), NGC 5194 (Sandage~ |l987| ), NGC 6946 (Freed- 
man et al. |1994|) and IC 342 (Krismer et al. |1995b. In or- 



der to take proper account of the mid-IR contribution to the 
FIR luminosities, the FIR luminosities were calculated from 
FIR = 5.35 10 5 r 2 (12.66 fi 2 +5.00 f 25 +2.55 f 60 +1.01 fi 00 ) L Q , 
where r is the distance to the galaxies in Mpc and fi 2 , f 2 5, 
f6o and fioo are the apparent flux densities in Jy (Belfort et 
al. 1987 1, instead of the normal approximation, e.g. FIR = 
3.75KF r 2 (2.58 f 60 +1.0 fi 00 ) L Q (Lonsdale & Helou |1985|). 



This choice makes a small but noticeable change in the distri- 
bution of the data in the diagrams. 

Fig 2 displays the distribution of Haro 1 1 and the compar- 
ison sample in the IRAS f2s/f60 vs. f2s/fioo diagram, show- 
ing the extreme position of Haro 1 1 . Fig. 3 shows the FIR/B- 
60/100 diagram for the sample galaxies. In the diagram are 
also indicated the loci of normal spiral galaxies and blue com- 



pact galaxies (Dultzin et al. 1988). As we see, Haro 11 also in 



this case lies at the extreme of the BCG distribution. The ex- 
tremely strong 25/i emission confirms the lack of cold dust and 
indicates the presen ce of a strong nuclea r starburst (Dultzin- 
Hacyan et al. 1988; Hawarden et al 1986; see also Boulanger 



Fig. 4 and 5 correspond to Malhotras et al. diagram 1 . To 
extract the most important global parameters involved in the 
relation reflected in fig. 4, we made a cross correlation test be- 
tween the [CII]/FIR flux ratio, the IRAS data (i.e. the ratios 
between the 12, 25, 60 and 100// fluxes), the FIR/B flux ra- 
tios and the FIR flux. We found strong correlations (> 6a) be- 
tween only three of these parameters: fgo/fioo, FIR/B and the 
[CII]/FIR ratio. 

In figure 4 the trend found by Malhotra et al. is confirmed 
and appears even stronger. But it is also evident that the posi- 
tion of Haro 1 1 deviates from the envelope defined by the other 
galaxies. The possible explanations for the deviant position of 
Haro 1 1 are either that there is an excess in the [CII] over FIR 
luminosity or that the f6o/fioo temperature is high, or both. Let 
us look at these two options one at a time. A consequence of 
lower metallicity is that the UV photons penetrate deeper into 
the cloud, thus increasing the volume of the PDR region at the 
expense of the colder molecular core. Strong observational sup- 
port of this in terms of a high [C II]/CO flux ratio have b een re - 
ported by e.g. Mochizuki et el. ( |1994| ), Poglitsch et al. ( |1995| ), 
Madden et al. (|1997[), Smith & Madden (|1997|) and Israel et al. 



(199f). Haro 11 behaves in the same way with its extremely 
high [CII]/CO ratio (> 4 10 5 ). The interesting question is how 
much the [CII]/FIR ratio will be influenced by varying metal- 
licity under these conditions. The models by Kaufman et al. 
show that [CII]/FIR is almost unaffected by metallicity varia- 
tions in cases where Go/n have normal values. This is also in 
agreement with observations. In 30 Dor and the cases discussed 
by Israel et al. however, the Go/n ratio is low (< 0.1) which re- 
sults in a diluted radiation field, an increased heating efficiency 
and brighter [CII] emission. In Haro 1 1 the Go/n ratio is nor- 
mal (Go/n=l cm 3 ), as seen from table 3. So, both observations 
and modelling argue against an increase in [CII]/FIR in Haro 
1 1 due to low metallicity. There is also another reason. In fig. 
5 Haro 1 1 adheres to the other galaxies in the diagram. There- 
fore as a first guess it is reasonable to assume that [CII]/FIR is 
normal while f6o/fioo is deviating. 

We concluded in section 3.2 that Haro 1 1 contains very lit- 
tle warm or cold dust. Haro 1 1 is a merging system of galaxies 
and it is possible that the merging galaxies have lost their cold 
dust components, e.g. extended disks, in the merger process. 
Is it possible that the lack of a 'normal' cold dust component 
could explain the high f6o/fioo ratio? Would the contribution 
from such a component be sufficient to move the position of 
Haro 1 1 into the region of the other galaxies? It is an impor- 
tant question to ask since it separates global phenomena from 
local ones. We have made a rough check of this possibility by 
selecting the galaxies in the comparison sample having appar- 
ent sizes sufficiently large to allow us to obtain IRAS data both 
from the very centre and also the integrated IRAS fluxes of 
the entire galaxy. In this way we can have an idea about how 
much the extended component influences the 60/100 ratio. It 
turns out that in the majority of the cases the IRAS temperature 
actually increases when the extended component is included. 
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Moreover, for the hottest cases, the hot component is always 
significantly stronger than the warm/cold component, if exist- 
ing. Our impression is therefore that the IRAS temperature of 
the warm dust of the most active regions in Haro 1 1 and the 
galaxies below Haro 1 1 in fig. 4 are similar. 

Although we expect a large scatter in fig. 4 we should ex- 
amine if there may be other parameters involved ruling the bal- 
ance between [CII] and FIR emission and explaining part of 
the scatter in the diagram and the deviant position of Haro 1 1 . 
Since both the spectrum of the radiation field and the spatial 
distribution of gas and dust in general differs between low- 
mass starbursts and galaxies with AGNs, we could suspect that 
we would also find systematic differences in the [CII]/FIR ra- 
tios between these types of galaxies. When comparing galax- 
ies classified in the optical region as Syl, Sy2, liners and HII 
galaxies we find that the HII galaxies appear to have a broad 
distribution in f6o/fioo while the Sy2 and liners tend to concen- 
trate towards high IRAS temperatures and high FIR/B. Simi- 
larly, in a statistical investigation of the IR/optical properties of 
BCGs as related to Seyfert and Liner galaxies (Dultzin-Hacyan 
et al. 1988, 1990 > it was shown that the f25/fioo is a powerful 



CD 



o 




index that can be used to discriminate between pure starbursts 
and cases where AGNs are important for the radiation field. 



-.8 -.6 -.4 -.2 .2 
log(f(60)/f(100)) 



Fig. 3. The FIR/B luminosity ratio vs. the IRAS temperature 
fgo/fioo for the same sample as in fig. 2. The full drawn box 
indicates the locus of normal spiral galaxies and the hatched 
bo x the l ocus of blue compact galaxies according to Dultzin et 
al. 
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3.6. What rules the [CIIJ/FIR ratio? 

Several possible physical mechanisms have been proposed to 
explain the trends in [CII]/FIR seen in fig. 4 and 5 (e.g. Malho- 



tra et al. 1997; Luhman et al. 1998). As we already mentioned 



Malhotra et al. argued that the reduced efficiency in photoejec- 
tion due to a hard radiation field could be a plausible explana- 
tion. The data for Haro 1 1 however, acts as a counterexample 
of this. Both optical spectroscopy and the 'hot' IRAS indices 
indicate a hard radiation field. Still, the [CII]/FIR is close to the 
values of normal, more passive galaxies. An alternative expla- 
nation, that have been discussed by Malhotra et al. and others is 
that [CII] actually is optically thick in many of the cases where 
the [CII]/FIR ratio is low. The general opinion seems to be that 
the [CII] line is optically thin. These conclusions are based on 
estimates of the optica l depth in a few well known cases, e.g. 
M82 (Crawford |1985| ). But M82 is behaving well in the dia- 
grams (see fig. 5) and the central column density of the gas is 
not abnormal. From the discussion above we see that, depend- 
ing on the velocity field and the relative amount of molecular 
gas in the C+ zone, the medium could become optically thick 
in [CII] at a column density of N ~ iq 22-23 cm -2 . Such high 
numbers can be found in the central regions of some luminous 



starburst galaxies. Gerin & Phillips (1998) elaborated on the 



low [CII]/FIR ratio in Arp 220 and it is interesting to compare 
this galaxy with M82 and Haro 1 1 (see fig. 5). They came to 
the conclusion that with an estimated column density of N ~ 
10 24 cm -2 , the dust opacity and emission would be sufficiently 
high to explain at least part of the low value of [CII]/FIR. But 
in fact there may be three sources responsible for reducing the 
ratio: dust opacity, dust emission and self absorption. What is 
maybe most important is how to explain why the [CII]/Lb ; ra- 



tio decreases with FIR/B. Therefore the combined effect of dust 
absorption and [CII] selfabsorption should be investigated. 

Fig. 6 shows the same data as in fig. 5, divided into different 
galaxy types, Syl, Sy2, HII and Liners. As we can see, there 
is a predominance of Sy2 galaxies and Liners at high FIR/B. 
Most of these galaxies are ULIGs. There are strong indications 
that Liners are galaxies with massive nuclear starbursts and 
likewise some of the Sy2 galaxies. Many of these galaxies are 
probably remnants of massive merges. Even though the lumi- 
nosity of some of them are not extreme, high column densities 
could possibly be obtained also in less massive nuclear star- 
bursts, provided they are sufficiently compact. Such very com- 
pact nuclear starbursts are preferentially found in the type of 
galax y we find at high FIR/B in the diagram (Heisler & Vader 



1986). In a simplified form fig. 6 shows the effect of increasing 



optical depth on the sample. Where the gas is optically thick in 
[CII] the ratio between the [CII]/B flux ratio should be approxi- 
mately constant. Where it is thin, [CII]/FIR should be constant. 

4. [CII] as a probe of high redshift starbursts 

As we can see from figure 4 and 5, the [CII] line in starburst 
galaxies has a maximum luminosity of a few per mille of the 
bolometric luminosity and constitutes the peak in the spectral 
distribution from radio to X-ray. It has therefore been proposed 
(e.g. Petrosian et al. |1969fc Loeb |1993|; Stark |1997[) that it can 



be used as a probe in searches for high redshift starburst galax- 
ies. At redshifts of z> 5-6, where current models of galaxy 
formation predict the first massive galaxies to form and where 
optical-near IR observations start to reach the confusion limit, 
the [CII] line is shifted into the submillimeter (submm) region 
and could be detected with a suitable submm telescope at a cold 
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Fig. 4. The [CII]/FIR ratio vs. the IRAS temperature f 60 /fioo 
for the same sample as in fig. 2. Upper limits are marked with 
arrows. 
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Fig. 6. The [CII]/FIR ratio vs. the FIR/B luminosity ratio as in 
fig. 5. Upper limits are marked with arrows. The different activ- 
ity classes have been coded according to: Circle-Syl, Square- 
Sy2, Triangle: Liner, Star: HII, Pentagon: Unclassified. The full 
drawn line shows an idealized approximation of the effect ot 
increasing optical depth with FIR/B. 

a present day L* galaxy at z=4-5 in a few hours, using a 10m 
submm telescope at the South Pole. At the formation epoch a 
massive galaxy would be several magnitudes brighter than to- 
day, i.e. ~ 10 12 Lq. Thus we would need to merge about 10 
Haro 1 1 type galaxies to produce this luminosity. A beam size 
of 20"-30" for a 10m telescope corresponds to about 100 kpc at 
the redshift we discuss. Thus it would easily encompass a com- 
pact merging group of dwarf galaxies and make a detection fea- 
sible. Moreover, with the advent of planned large sub-mm/mm 
arrays (e.g. ALMA), the prospects for using the [CII]158 line 
as a probe of high redshift starburst would be even more inter- 
esting. 



Fig. 5. The [CII]/FIR vs. the FIR/B luminosity ratios for the 
same sample as in fig. 2. M82 is marked with a triangle and 
Arp220 with a square. Upper limits are marked with arrows. 



dry site like the South Pole (Stark 1997). The most interesting 
targets would probably be massive starburst mergers, character- 
ized by a high 60/100 ratio and a high FIR/B ratio. The tenden- 



cies seen in the study of Malhotra et al. (1997) however, have 
created some pessimism (e.g. Gerin & Phillips 1998} ) about the 
feasibility, due to the rapid decrease in the relative strength of 
the [CII] line with IRAS temperature. The properties of Haro 
1 1 however, indicate that the conditions may be more promis- 
ing than anticipated since the [CU]/Lb i of this galaxy is close 
to the value obtained for the normal galaxies. If this is due to 
the low metallicity of the galaxy or not remains to be explored 
but if so it may apply on distant met al-poo r galaxies. The detec- 
tion probabilities derived by Stark ( 1997 ) predict that it would 



be quite possible to detect a galaxy with a mass equal to that of 



5. Conclusions 

We have reported about ISO LWS observations of the luminous 
metal poor BCG Haro 1 1 in the three prominent emission lines 
[OI]63/i, [OIII]88^ and [CII] 1588//. The galaxy is involved in 
an intense global starburst probably due to a merger of dwarf 
galaxies. The ISO data and previous observations in the optical- 
radio region in combination with model predictions are used to 
derive information about the physical properties of gas and dust 
in the galaxy. 

Haro 1 1 is one of the hottest IRAS galaxies observed and 
we find no trace of cold dust. Most of the neutral hydrogen 
is confined to the photodissociation regions. The main part of 
the gas however, appears to be in ionized or molecular state. 
This possibility is generally neglected when discussing the gas 
content in starburst dwarfs, leading to serious underestimates 
of the total gas mass. 

We also reinvestigate the claimed correlation between [CII] 
emission, far-IR luminosity and the IRAS 60/100 temperature 
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for a number of starforming galaxies of different types, active 
and non-active. We confirm and tighten the correlation and note 
that Haro 1 1 deviates in the sense that its [CU]/Lfir, is higher 
than expected. We argue that this may indicate that the pre- 
viously preferred explanation of the relationship, a decreasing 
efficiency with IRAS temperature in the ejection of photoelec- 
trons from UV-illuminated grains may not be entirely correct. 
An alternative, or complementary explanation seems to be that 
the optical depth increases with increasing IRAS temperature. 
This would be expected if the IRAS temperature correlates with 
compactness and/or mass of the central starburst region. The 
metallicity may also play an important role, in which case the 
somewhat pessimistic opinions about using [CII] 158/i to detect 
high redshift massive starburst galaxies may not be valid. The 
whole situation will probably be clarified once a more repre- 
sentative sample of galaxies, including both normal dwarfs and 
massive galaxies becomes available. 
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